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Abstract 

The bias dependent interface charge is considered as the origin of the ob- 
served non-ideality in current-voltage and capacitance- voltage characteristics. 
Using the simplified model for the interface electronic structure based on de- 
fects interacting with the continuum of interface states, the microscopic origin 
of empirical parameters describing the bias dependent interface charge func- 
tion is investigated. The results show that in non-ideal metal-semiconductor 
contacts the interface charge function depends on the interface disorder pa- 
rameter, density of defects, barrier pinning parameter and the effective gap 
center. The theoretical predictions are tested against several sets of published 
experimental data on bias dependent ideality factor and excess capacitance in 
various metal-semicoductor systems. 



1 Introduction 

The Schottky barrier heights at the metal-semiconductor interfaces are usually found 
to be only weakly dependent on the metal work function pQ. The failure to comply 
with the Schottky-Mott rule jSIHj, where the Schottky barrier height, is equal 
to the difference between the metal work function, m , and semiconductor electron 
affinity, Xs, is often termed Fermi level pinning. This phenomenon has stimulated 
extensive debate over the past few decades, resulting in several proposed models, and 
has recently again attracted considerable attention in the scientific community [3J 
IZj- Interface properties enter these models in the form of the charge localized 
near the interface resulting from the finite density of states in the band gap and/or 
defect levels. It has been recently shown |H1 El Ej that delibarate introduction of 
metal atoms into the thin semiconductor layer near the metal-semicoductor interface 
utilizing ionized cluster beam deposition [TUl HU C2 E] , enables identification of the 
origin and analysis of non-ideality in certain metal-semiconductor contacts. 

The non-ideality is a consequence of the interlayer present in the structure between 
the metal and the semiconductor substrate. In ionized cluster beam deposited sam- 
ples, this interlayer consists of the penetrated metal atoms embedded in the host 
semiconductor lattice where they introduce disorder into the lattice and may act as 
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additional deep impurities. The incident beam in ionized cluster beam deposition is 
comprised of neutral and, to smaller extent, also ionized metal atoms [TTH ITT1 IT2] . 
The small number of neutral and ionized metal clusters of various sizes which are si- 
multaneously present in the incident beam, on account of the very low energy shared 
by each cluster atom (of the order of about 60 eV per atom for an avarage cluster 
size of 16 atoms and 1 keV translational energy) can be safely neglected. The addi- 
tionally gained energy, up to few keV in some cases, enables the metal ions to enter 
into semiconductor bulk and come to rest within and up to their penetration length 
(L) typically a few nm in depth. In such a non-ideal metal-semiconductor struc- 
ture one can identify two distinct interfaces, one situated at the metal/disordered 
interlayer, and the other at the disordered layer/regular semiconductor junction 
It was shown that the electronic properties of such a non-ideal Schottky structure 
determine its low-frequency capacitance and Fermi level pinning behaviour jHllHl- 

Closely connected to the questions of the microscopic properties of the m/s interface 
is the problem of the excess capacitance O EH EH EH El HI in Schottky junctions. 
Recently [S] we have been able to show that the observed excess capacitances in 
non-ideal Schottky junctions arise on account of the bias dependent localized charge 
density induced at the appropriate interface in the structure. The model function 
for the interface charge was chosen to be: 

a(U) = bexp[-(K + qU/U ) 2 ], (1) 

where the constants b, K, and Uq were determined by comparison with the exper- 
imental capacitance- voltage dependence of different systems |S]. The differential 
capacitance is in this model given by 

where L is the width of the disordered layer, and W is the width of the depletion 
layer. Here, the dielectric constants of the disordered layer and the homogenuous 
semiconductor (e s ) are taken to be equal. As an example, we show on figure 1 the 
measured low frequency capacitance (plotted as C~ 2 {U)) of e-beam deposited Al/n- 
Si (data taken from refs. [d EB]) and of ICB deposited Pb/p-Si (our own samples) 
Schottky structures as a function of applied bias. The full line is the calculated 
capacitance-voltage characterstics based on eqs. E] and El For further details we 
refer the reader to reference 

The current-voltage characteristics of non-ideal metal-semiconductor contacts do 
not usually follow the characteristics predicted by thermionic emission transport 
theory. The observed discrepance is often treated by introducing the ideality factor 
into the current- voltage characteristic defined by [2|: 

where (ftb is the barrier height. 



In what follows we show that the bias dependence of the ideality factor used to 
describe the non-ideality of the current-voltage characterstics can be traced back to 
the interface charge causing the non-ideality and later on we discuss the microscopic 
origin of the parameters in the interface charge functions. The derived expressions 
are used on experimental sets of data published by different authors, and some of 
the values of microscopic parameters are extracted. 



2 Ideality factor and interface charge in non-ideal 
Schottky contacts 

In a non-ideal metal-semiconductor system we can quite generally describe the non- 
ideality with the macroscopic charge situated near the metal-semiconductor inter- 
face. The charge neutrality of the system demands: 

Qm + Qsc + Qi = 0, (4) 

where Q m , Q sc , and Qi are the surface charge densities on the metal in the depletion 
region of the semiconductor, and at the interface. 

For the non-ideal metal-semiconductor structures with the disordered interlayer and 
the bias dependent charge density the charge neutrality is written as: 

Q m + Qsc + Qi = Qm + qn 2 W + q{n 2 + n x )L + a{U) = 0, (5) 

where n 2 is the concentrations of the ionized bulk semiconductor dopants, n\ is the 
concentrations of the ionized defects in the disordered interlayer, and L is the width 
of the disordered interlayer. The bias dependent localized charge density residing 
at the disordered layer /homogenous semiconductor interface is here introduced with 
Pl = a(U)5(z-L). 

The electrostatic analysis of the model yielded the following expression for the width 
of the depletion layer: 

W = Jtf (-L -qU + fe . o - Xs) , (6) 

where U is the applied bias, 4>b,o is the Schottky barrier height (potential barrier at 
the metal/disordered layer interface). 

The potential at the disordered layer /semiconductor interface is defined as the ef- 
fective Schottky barrier height in non-ideal contact: 

<P b = ^(L-W) 2 + qU + Xs- (7) 

Using this and the definition of the ideality factor (eq. ([3}) we arrive at the following 
expression for the bias dependent ideality factor in non-ideal metal-semiconductor 
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systems: 



1 W(U)J V e s dU 



This result limits to the well known result of Card and Rhoderick ^3] when L « 
W . If we introduce depletion and interlayer geometric capacitances Cd = e s /W ', 
Ci = e a /L, we obtain: 

<U) = (»> 

1 "t" ddU 

Card and Rhoderick [THj proposed the division of the interface states into two groups: 
one in equilibrium with the metal, and the other in equilibrium with the semiconduc- 
tor. Gomila and Rubi UHj have considered the problem of the bias dependence 
of ideality factor in case of nonequilibrium, and using a kinetic model derived ap- 
propriate bias dependence of these partial densities of states. 

The detailed electronic structure of the non-ideal metal-semiconductor contact must 
be complex involving possible defect, metal and disorder induced densities of states. 
The charging mechanisms at the interface are therefore quite complicated. For 
instance, Wu, Yang, and Evans [21] described the observed excess and negative ca- 
pacitance in NiSi 2 /n-Si(100) Schottky diode using the modified occupation function 
of the interface states. They introduced modification causing the decreasing of the 
interface charge with the forward bias on account of the decreasing of the occupancy 
of the interface states below the metal Fermi energy. 

Here we use a much simpler model for the complex electronic structure which how- 
ever gives suprisingly good description of the functional bias dependence of the 
excess charge. The two defect levels are considered as two sharp electronic levels e\ 
and €2 interacting with the partial densities of states in equilibrium with the metal 
and semiconductor. This mechanism is in some respect similar to the description of 
the interaction of the metal clusters on the semiconductor surface [22] • The impor- 
tant difference is that in our case we consider a single defect coupled to the interface 
states, so that the transferred charge cannot be screened out as it is the case in the 
metal cluster [22]- 

The random distribution of the metal atoms in the disordered interlayer causes 
the fluctuations of the potential at leads to broadening of the sharp defect levels. 
In the high density limit the distribution function of the potential V(r) can be 
approximated by the gaussian: 

P{V) = -L-eM-V 2 /(U )% (10) 
V nU o 

where the parameter Uq from the charge model function was immediately recog- 
nized as the width of the distribution function. The disordered density of states in 
semiclassical approximation yields [~ 



Pi(e) = ^-(exp(-(e- ei )/^o) 2 )). (11) 
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P2(6) = - 7 ^-(exp(-( e -e 2 )/f/o) 2 )). (12) 

If the impurity density is n then we obtain for the excess charge density using the 
low temperature limit (T — ► 0) for simplicity: 

2/3 

a{e f ) = ^- (erf [(e f - (e d - Ae))/U ] + erf [(e d + Ae - e f )/U \) , (13) 

where e/ is the semiconductor Fermi energy, n is the defect concentration, is the 
mean defect level, and 2Ae is the separation between the two defect levels. If the 
shift of the defect levels from the mean value is small compared to the disorder 
parameter Uo, i.e. Ae/Uo « 1, then we have approximately: 

2*777 2 / 3 Af 

« " 1 7^exp[-(6 / /C/o - e d /U f]. (14) 

The difference between the Fermi energy in the metal and the Fermi energy in the 
semiconductor far from interface equals the applied bias: 

e f = e fm + qU. (15) 

In equilibrium, the Fermi energy at the interface is found to be pinned near the 
midgap of the semiconductor. In one dimension it tends to fall near the branch 
point in the complex band structure. In three dimensions the pinning point is found 
to be given by the gap center of the minimum indirect gap, e g , corrected with spin- 
orbit splitting [25J, and including the shift from the defined center depending on the 
metal 5 m : 

e/ m = -e 9 + £ m . (16) 

Using (4), and (5) in the expression for the interface charge density (3), one can 
define an effective gap center, eo,e//, in the presence of the disorder and defects in 
the semiconductor: 

1 1 

e 0>e ff = KU = -e g + 5 m - e d = -e g - 5e m>d . (17) 

This shows that the parameter K of the model charge function (1) identifies the 
position of the effective gap center at non-ideal metal-semiconductor interface. The 
Schottky barrier height is in most models describing the Fermi level pinning 510311] 
given as an interpolation between the Schottky-Mott and Bardeen limit: 

<fb, P = c(I s - (j) m ) + (1 - c)e 0) eff. (18) 

Here, I s = E g + \s is the semiconductor ionization potential, and c is the pinning 
parameter which depends on the details of the pinning mechanism. As the pinning 
parameter changes its value from 1 to 0, the Schottky barrier height goes from totally 
unpinned to fully pinned value. 
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3 Results and discussion 



Using the experimentaly observed values for the Schottky barrier height on p-type Si 
and the effective gap center deduced from the analyzed capacitance-voltage charac- 
teristics jHj one can obtain the pinning parameter for different metal-semiconductor 
systems that were analyzed If the pinning mechanisms depends on the contin- 
uum of the localized density of states in the band gap, then the pinning parameter 
of the following form is usually obtained |2j for nearly ideal systems: 

(19) 



j _|_ q 2 D(e , e ff 



where _D(eo, e //) is the (two dimensional) density of states in the band gap at the 
effective gap center, £ is the localization length of the states, and e s is the dielectric 
constant of the semiconductor. 

When the disordered interlayer of width L is introduced between the metal and 
the bulk semiconductor, possibly extending over several atomic layer, the pinning 
parameter is given by 0: 



c = sech \ Jq 2 L 2 p(e ,eff)/e s ) , (20) 
where p(eo,e//) is the density of disorder induced gap states. 

The energy interval determining the properties of the non-ideal metal-semiconductor 
junction is estimated to be of the order of Ae. The density of the defect levels 
interacting with the continuum of localized states is then approximately given by: 

D(e , eff )Ae « n^. (21) 

Using (8) and (9) in eq. (3) gives the final expression for the interface charge density 
in nearly ideal systems: 

" (U) = (l/c-l)Je,U» eM - (e °-" /U ° + " U/U ° n <22) 
and the corresponding expression for the non-ideal systems: 

= (m cJ(c)) W . «Pt-( + tfWfi (23) 

The expression for the interface charge density (eq. 10) as derived from the pre- 
sented model shows that the crucial microscopic parameters detemining the interface 
charge density in non-ideal metal- semiconductor contacts are: the disorder param- 
eter Uq, density of defects n, Schottky barrier pinning parameter c (or the density 
of states interacting with the defect level and localization length), and the effective 
gap center. The parameters b (the maximum of the interface charge density) and 



6 



Uq obtained from the analysis of the capacitance- voltage characteristics of non-ideal 
metal-semiconductor contacts jH] allows one to compute the number of defects per 
unit area, the density of corresponding continuum of states and the shift of the 
defect level interacting with the continuum of localized states. The results of the 
calculation are presented in tables 1 and 2. The calculated effective gap center for 
the series of metals on p-type Si (Ag, Pb, Al, Ti, Mo) reasonably well corresponds 
to the calculated avarage branch point, e = 0, 36 eV, in the Si complex band struc- 
ture > the maximum deviation being 0,1 eV. 

We have also considered the negative and excess capacitance observed by Wu, Yang, 
and Evans |2H i n their samples. The excess capacitance is in the case L « W 
according to eq. El simply equal to derivative of the interface charge with respect to 
applied bias: C excess ~ —da/dU. Using generalized model for the bias dependence 
of the interface charge eq. ( IT^ : 



a(U) = gmerf (e + Ae - qU)/U 0>1 ) + qn 2 erf (qU - e + Ae)/tf , 2 ) , (24) 



we are able also to describe the results of Wu, Yang and Evans [21] within the 
framework of our model as it is shown in figure 2. The parameters used in calculation 
were: e = 0, 5 eV, Ae = -0, 1 eV, U 0A = 0, 065 eV, U 0j2 = 0, 17 eV, n x = 3, 8 x 10 14 
cm" 2 and n 2 = 3, 1 x 10 14 cm -2 . 

We now apply the derived results for the bias dependent ideality factor given by 
eq. ( IHJ) to two sets of published data measured in quite different non-ideal metal- 
semiconductor systems. 

Maeda et al. [2E] have observed bias dependent ideality factor in Au/n-GaAs metal- 
semiconductor system. The bias dependence of the ideality factor exhibited a pro- 
nounce peak at certain applied bias. We have used eqs. ( HJ) and ( |HI) and calculated 
the measured bias dependence of the ideality factor using the following set of pa- 
rameters: b = 4,6 x 10~ 3 As/m 2 , U = 0,33 eV, K = 0,042, L = 2 nm. Figure 3 
shows the measured set of data and the calculated bias dependence of the ideality 
factor. In tables 1 and 2 we present the computed density of defects, the shift of 
defects energy levels and the density of states. 

In figure 4 we present the experimentally observed bias dependence of ideality factor 
as measured in Au/n-InP Schottky structure obtained from [21| (empty dots). In 
the applied voltage interval for the presented set of data (2Zj we can write for the 
derivative of the charge density: 



with the following values of parameters: c\ = 0, 5, c 2 = Uq = 0, 18 eV, C3 = e^ = 0, 11 
eV, c 4 = 0,25 1/V, c 5 = 0,19. The width of the interlayer was L — 15 nm |27| . 
From these values we calculate the density of the defects = 2, 8 x 10 11 cm -2 and 
the barrier height for Au/InP contact: ^ „ = 0, 76 eV which agrees well with the 
results obtained from analysis of I-V curves in ref. [2*7j . 




(25) 
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4 Conclusions 



In this work we have shown that the non-ideality in metal-semiconductor contacts 
which is observed either in current-voltage as bias dependent ideality factor or in 
capacitance voltage characteristics as excess capacitance can be described using a 
single model for the bias dependent interface charge. 

The expression for the interface charge density is derived which showed that the 
crucial microscopic parameters detemining the interface charge density in non-ideal 
metal-semiconductor contacts are: the disorder parameter U , density of defects n, 
Schottky barrier pinning parameter c (or the density of states interacting with the 
defect level and localization length), and the effective gap center eo, e //- The param- 
eters b (the maximum of the interface charge density) and Uq obtained from the 
analysis of the capacitance- voltage characteristics of non-ideal metal-semiconductor 
contacts jH] allows one to compute the number of defects per unit area and the 
shift of the defect level interacting with the continuum of localized states. The 
basic assumption was to condsider the random distribution of the metal atoms in 
the disordered layer causes the fluctuations of the potential in the layer resulting in 
broadening of the sharp defect levels. In the high density limit where the distribution 
function of the potential can be approximated by the gaussian. 

The bias dependent charge residing at the metal- semiconductor interface is con- 
nected to the bias dependence of the ideality factor of the cur rent- volt age character- 
istic. The comparison of the results computed from the theoretical model with the 
published experimental data enabled us to extract several microscopic properties of 
various metal-semiconductor contacts. 
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Figure 1: Low frequency capacitance of e-beam deposited, Al/n-Si as a function of 
the external reverse bias [T71 [TH] , and low frequency capacitance of ICB deposited 
Pb/p-Si(100) as a function of the external reverse bias are compared to the calculated 
values jHJ. 
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Figure 2: Values of ideality factor (empty dots) in Au/n-GaAs(100) as obtained 
from experiment j2E], and calculated bias dependence of ideality factor (solid line). 
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Figure 3: Measured low frequency susceptance as a function of applied bias (empty 
dots) of e-beam deposited NiSi2/n-Si(lll) [13121] and calculated susceptance (solid 
line) using the two-level model for the bias dependent interface charge. 
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Figure 4: Values of ideality factor (empty dots) in Au/n-InP(100) as obtained from 
experiment j22j, and calculated (solid line) bias dependence of ideality factor. 
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Table 1: Values of experimentaly determined effective gap center, barrier height and 
calculated values of the pinning parameter and the density of states. The values for 
the second column were taken from J2[ E^, and for the third and fourth column 
from L 8l. In the calculation the localization length was taken to be 0,25 nm for 
the nearly ideal metal-semiconductor contacts. In non-ideal cases the width of the 
interlayer was used. 
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Table 2: Values of experimentaly determined maximum of the interface charge den- 
sity b, disorder parameter, and calculated defect density and energy level shift. 
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